Electronic Transport in Single-Molecule Magnets on Metallic Surfaces 
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An electron transport is studied in the system which consists of scanning tunnehng microscopy- 
single molecule magnet-metal. Due to quantum tunneling of magnetization in single-molecule mag- 
net, linear response conductance exhibits stepwise behavior with increasing longitudinal field and 
each step is maximized at a certain value of field sweeping speed. The conductance at each step 
oscillates as a function of the additional transverse magnetic field along the hard axis. Rigorous 
theory is presented that combines the exchange model with the Landau-Zener model. 

PACS numbers: 75.45.-|-j, 75.50.Xx, 75.50.Tt 
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Recently high-spin molecular nanomagnets such as 
Mni2 or Feg attracted lots of attention due to obser- 
vation of quantum tunneling of magnetization and pos- 
sible applications in information storage and quantum 
computing 1, 2, 3, 4, 5, 6]. These single-molecule mag- 
nets (SMMs) exhibit steps in the hysteresis loops at low 
temperature, which is attributed to resonant tunneling 
between degenerate quantum states or quantum tunnel- 
ing of magnetization(QTM). These unique features of 
SMMs are the consequence of long-living metastable spin 
states due to the large spin and strong anisotropy of 
SMMs. QTM also made it possible to detect the inter- 
ference effect of Berry's phase on the magnetization at 
each step while the transverse field along the hard axis 
is variedj^l^. Novel features of quantum tunneling are 
expected to manifest themselves in, if any, other observ- 
ables. Especially the effects of QTM on the electronk 
transport remain to be explored in both experiments 7] 
and theories. 

In this paper we study theoretically the effects of QTM 
on the transport properties of SMMs which are deposited 
on a metallic surface with monolayer coverage. Placing 
the scanning tunneling microscopy(STM) tip right above 
one SMM, we compute the electric current which flows 
through a SMM when the bias voltage is applied between 
the STM tip and the metallic substrate(Fig. P). We find 
that the linear response conductance increases stepwise 
like the magnetization of a SMM as a longitudinal mag- 
netic field is increased. The stepwise behavior of conduc- 
tance results from the QTM in SMM. The conductance 
at each step oscillates periodically as a function of addi- 
tional transverse magnetic field along the hard axis. Our 
theoretical predictions are not known in the literature as 
far as we know and can be tested experimentally. 

When a finite bias voltage is applied between the STM 
tip and the metallic substrate, the electrons will tunnel 
through a vacuum between the metal surface and the 
STM tip. Since the STM tip is placed right above the 
SMM in our model system, the tunneling electrons may 
well be scattered by the large spin of a SMM. Our model 
system can be considered as the conventional tunnel junc- 



tion with a SMM sandwiched between two normal metal- 
lic electrodes. The metallic substrate and STM tip are 
conveniently called the left(p = L) and right (p = R) 
electrodes, respectively. Two electrodes are described by 
the featureless conduction bands with the energy disper- 
sion epk, Hp = J2ka '^pk<^\ka%ka- ^hc Hamiltouian of the 
SMM will be introduced later. Tunneling electrons are 
modeled by the HamiltonianP 



kk' a. 

+ (;^LBAkJ<^P'^Rk'p ■ -S' + H.c.) , (1) 

ka k' p 

where a and f3 indicate the spin direction of electrons. 
The first line represents the direct tunneling between two 
electrodes, while the second line describes the tunneling 
of electrons scattered by the spin S of SMM. Our the- 
ory is equally applicable to the molecular break junction 
geometry. 

The electric current can be computed using the 
Keldysh Green's function method or equivalently the 
Fermi's golden ruleQ. In this paper we study the very 
weak coupling limit so that the higher order process like 
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FIG. 1: Schematic diagram of our model system. A single- 
molecule magnet (SMM) is deposited on a metallic surface 
and the scanning tunneling microscopy (STM) tip is posi- 
tioned right above the SMM. The easy axis of SMMs is di- 
rected normal to the metallic substrate. 
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the Kondo effect may be safely neglected. In this case it 
is enough to compute the electric current up to the lead- 
ing order term. Using the Fermi's golden rule the electric 
current can be written as 

Ilr = P„i ^ ^ VKLfeam^_R.fe'/3m' 

m ka k'l3 

x/(eL/c)[l - fi^Rk')] - (Lkam ^ Rk' (3m'){2) 

Here Wi^j is the transition rate from the state i to j, 
/(e) is the Fermi-Dirac distribution function and Pm is 
the probability for the SMM to be in the state Sz = m. 
The leading contribution to the transition rate is given 
by the expression W,^j = ^| < > \^S{E, - E,), 

where i and j are the collective indices denoting the states 
{Lkam} or {Rk' (3m'} and Epkam = Spk + Em with epk = 
{p = L,R). Hp is the chemical potential shift in 
the electrode p due to the source-drain bias voltage, and 
Em is the energy of the state 5*2 — m in the SMM. 

Up to the second order in Tlr and Jlr, we find the 
electric current to be 
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where 7t(7j) = '^'^'^NlNb\Tlr\^{\Jlr\^) is a measure of 
the dimensionless direct (spin-scattered) tunneling rate, 
(•^z) = Sm"^^-^™' ^ source-drain bias voltage 

given by eV = hl - Mfl: and C(e) = e/[l - exp(-/3e)] 
with I3~^ — ksT. The linear response conductance is 
then G = 2|! [^^ + ^^g,(T)] , where g,(r) = (S^^) + 

J2^Pm[SiS+l)-m{m±lMEm-Em±i) with 77(e) = 
dC{e)/de. We would like to emphasize that only the spin- 
exchange tunneling reflects the dynamics of the QTM 
inside the SMM. 

Due to the crystal electric field arising from the struc- 
ture of a magnetic molecule, the ground state spin mul- 
tiplet cannot remain degenerate. The effective Hamilto- 
nian for the ground state spin multiplet of independent 
SMMs such as Fes can be expanded as|3.llo| 



■HsMM = -DSl + E{S^-S^y) + C{S[ 
-giig{HzSz + H^Sx), 



(4) 



where Sx, Sy, Sz are three components of the spin op- 
erator, S± = 5*3, ± iSy, D and E are the second-order 
and C the fourth-order anisotropy constants, and the 
last term is the Zeeman energy. In the absence of 
transverse terms, the energy level of the state Sz — m 
is Em — —Dvn? — g^gHzTn. When we start with a 
ground state Sz — —S corresponding to a large neg- 
ative longitudinal field, the level crossing with states 
Sz = S — M {M = 0, 1, 2, • • • ) occurs at resonant fields, 
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- MD/gfiB- When Hz 
states Sz = —S and Sz = S — M are degenerate energet- 
ically. Turning on the transverse terms leads to mixing 
of two degenerate states, lifts the degeneracy at the res- 
onant fields and results in the avoided level crossing. 

The scaled conductance gs can be simplified as 
gs{M) = S"^ + X^nLo ""^■S-n at zero temperature by 
noting that Es < Es-i < ... and 77(e) — 9{e), the 
step function. In deriving this expression of gs{M) it 
is assumed that the weight transfers from Sz = —S to 
Sz = S, S — 1, ■ ■ ■ with increasing longitudinal magnetic 
fields. 

To compute the probability, we need to solve the 
time-dependent Schrodinger equation for the Hamilto- 
nian Hsmm- The probability is defined as Pj = 
limt_^oo |fij(i)P when the wave function is written 
as \'9{t)) = J2j=-s "^^^ time-dependent 

Schrodinger equation for \^{t)) is reduced to the coupled 
25+1 differential equations for the coefficient aj{t). Re- 
cently it was numerically found [lollll| that the two-level 
approximation can reproduce quite well the results of the 
full differential equations. In the ensuing discussion we 
adopt the two-level approximation to find an analytic for- 
mula of the probability. The weight transfer is found to 
occur only between the states Sz = —S and Sz = S — M 
at the resonant field H^fj', for M = 0,1,2,.., until the 
complete depletion of the state Sz = —S. The amount 
of such weight transfer depends on the magnitude of the 
tunnel splitting or mixing Am between two states. At 
the resonant field H^j^ the full Hamiltonian H smm is 
approximated as the effective two-level model|lolll2| be- 
tween the states Sz = ~S and S^ ~ S — M, 
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where c(= dHz/dt) is the field sweeping speed. 

Defining r = gnBct/AM and 7^ = hg^iBc/A\;, we 
find the coefficient [l^ around the resonant field H^m ^ 
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where aM = {2S - M)/(27m), Am = l/(8aM7M): 
Fj — exp(— 27rAj) and D is the parabolic cylinder 
function[lJ|. The desired probabilities are then Ps-m = 

(1 - FMmfr^' Fj) and P_s = 11^=0 ^3- Note that F, 
and 1 — Fj denote the probability for an SMM not to 
transfer and to transfer from Sz~—S to Sz — S — jaX 
the j-th resonant field, respectively. 

To illustrate the above analytical results with con- 
crete number, we compute the scaled conductance, gs{= 
gs — S"^) at zero temperature for an octanuclear iron(II) 
oxo-hydroxo cluster of formula [Fe802(OII)i2(tacn)6]^~'' 
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FIG. 2: The scaled conductance gs{M) vs. the longitudinal 
field Hz at zero temperature for three typical sweeping speed 
(T/sec). M = 1, 2, 3, 4 indicate the positions of the resonant 
fields, H^i^^ 0.215, 0.429, 0.643, 0.858. in units of Tesla. 
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FIG. 3: Dependence of ga{M) on the field sweeping speed c 
at each resonant field. Inset: the magnetization (Sz) vs. c at 
each resonant field. 



where tacn is a macrocyclic ligandfl. We adopt the 
model parameters from Refs. |a and [lOf D = 0.292K, 
E = 0.046K, C = -3.2 x IQ-^K. The tunnel splitting 
Am is calculated for = O.lHz at the resonant field 
by employing the numerical diagonalization 10| or the 
perturbation methodT^. We obtain qualitatively the 
same results when has the fixed value at all resonant 
fields 0. 

The scaled conductance, gs{M) = X]f=i Ilj^o ~ 
A/n*lo ^hich is vahd for uf^ < < H^m+i, is dis- 
played in Fig. El for three typical field sweeping speeds. 
Similar to the magnetization curve, the scaled conduc- 
tance is featured with the stepwise increase as a function 
of magnetic fields. The jumps in gs{M) occur at the res- 
onant fields and are caused by the QTM in SMMs. The 
step height is very tiny (- 0.318 x 10"^) at uf^ = 0.215 
T for all three sweeping speeds. At the second and third 
resonant fields the step heights are more pronounced and 
their magnitude depends sensitively on the value of c. 
Some steps are missing depending on both the sweeping 
speed and the resonant fields. 

To study in more detail the structure of the steps 
in the conductance we plot in Fig. |31 the scaled con- 
ductance gs{M) at each resonant field as a function of 
the sweeping speed c. In comparison the magnetization, 

(^.> = S- Y.ti n:=o F, - {2S- M) n;io dis- 
played in the inset. The magnetization is a monotoni- 
cally decreasing function of c while the conductance is 
nonmonotonic and maximized at the specific value of c. 
Since the weight transfer, 1 — Fj at h'"^\ from Sz = —S 
to Sz = S — j is monotonically decreasing with increasing 
c, the magnetization is expected to decrease with c. 

Unlike the magnetization, the conductance has con- 
tributions only from the transferred states but not from 
Sz = —S. Since Fj is increasing with c, P_5 is an in- 
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FIG. 4: Oscillation of qm as a function of transverse field for 
c = 0.014T/sec. For this sweeping speed the M = 3, 4 curves 
are almost identical except — 0, 0.4, and 0.8 T. 

creasing function of c while Ps~m has the maximum 
value as a function of c. The conductance gs{M) has 
the contribution Sgs = MPs-m from the Af-th reso- 
nance and is expected to have the maximum value at 
some value of c. Such a sweeping speed can be com- 
puted approximately as c"^^^^^ ~ [i: / {2hg^B)][^\i / {"^S — 
Af)][log[MEf£o^^/(l + E*LT'ELo^O]]-^ where = 
(2S')Af/[(2S' - i)A2]. The values of c^-^""""'^ (T/sec) are 
5.1 X 10"^ 1.08 X 10"^ 0.182 at M = 1, 2, 3, respectively. 
Even though there exists a maximum in the scaled con- 
ductance at M — 4, the value of c = 5.16 (T/sec) lies 
beyond experimentally meaningful range. In order to 
observe the steps in conductance at Af = 3 or Af = 4 res- 
onance, the sweeping speed should be larger than about 
0.01 or 0.1 (T/sec), respectively. 

The conductances at the resonant fields are displayed 
in Fig. ^ as the transverse field is varied along the hard 
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axis. Similar to the magnetization the conductance at 
each resonant field oscillates with almost the same pe- 
riod of ~ 0.4 T. Such oscillatory conductance faithfully 
reflects the structure of the tunnel splittings as a func- 
tion of the transverse field^^. The periodic modulation 
of tunnel splittings by the transverse field results from the 
interference between two spin paths of opposite windings 
around the hard axis [111,113, 111 . 

The tunneling splitting is known to vanish at the lat- 
tice of the diabolic fields At such fields the tunneling 
probability is zero so that the jump in the conductance 
vanishes. Depending on the parity of M, the oscillations 
of the conductance have the different phase. The M = 2 
curve is out of phase compared to the M — 1,3 curves. 
For example, the conductance for M = 2 takes on the 
minimum value at the transverse field where the conduc- 
tance for M = 1 is maximized. This parity behavior 
originates from the impossibility of matching an even- 
valued wave function with an odd- valued one which gives 
rise to diabolic fields. Weak structures around — 0, 
0.4, and 0.8 T for M = 3,4 curves can be made con- 
spicuous with varying the field sweeping speed. Though 
the overall structure of oscillatory conductance persists, 
the amplitude of oscillations depends sensitively on the 
sweeping speed 16]. 

We briefly address the effect of experimentally relevant 
issues on our theoretical results. It may be important to 
consider the effect of environmental degrees of freedom 
such as phonons, nuclear spin and dipolar interaction 
[l^ on the magnetization process of SMMs. Such in- 
teractions make the SMM relax to the true ground state 
Sz = S and the relaxation process helps the magneti- 
zation to recover its full stretched value. Since all the 
transferred states Sz = S - M{M = 1, 2, • • • ) lose the 
weight to the ground state, we expect that the value of 
Qs will rise stepwise with increasing field and might van- 
ish in the end due to the relaxation process. Since the 
elapsed time between steps, which is of the order of 10 
sec or less for the typical sweeping speeds (see Fig. ^ , 
is much smaller than the relaxation time of magnetiza- 
tion (- 10* sec) H [13, we believe that the stepwise be- 
havior of the conductance can be observed experimen- 
tally in the typical field sweeping speeds. The effect of 
anisotropy in SMMs on the conductance was clarified in 
our work. In the absence of anisotropy gs = S{S + 1) so 
that the anisotropy in SMMs modifies the conductance 
by the amount S out of S{S -1-1). In the case of Fes 
or Mni2 5 = 10 so that the modified conductance is es- 
timated to about 10% which lies in the experimentally 
detectable range. Possible exchange anisotropy in spin- 
scattered tunneling can be addressed, 16j by considering 
the ratio, a = (J£^ + J^j^) V4[J|^]2.' When a > 1, the 
conductnce steps are more enhanced than the isotropic 
case(a — 1). For the case of a < 1, the steps are reduced 
or can be negative depending on the value of a. 

In summary we studied the current-voltage character- 



istics of the STM-SMM-metal system at low tempera- 
ture. We found that the quantum tunneling of magneti- 
zation (QTM) in SMMs has a substantial effect on the 
electronic transport. The QTM in SMMs leads to the 
stepwise behavior in the conductance (just like the mag- 
netization) when the magnetic field is applied along the 
easy axis. Unlike the magnetization the conductance at 
each resonance is nonmonotonic with the sweeping speed 
and reaches the maximum at some sweeping speed. In 
addition, the conductance at the resonant fields is oscil- 
lating as a function of the transverse field applied along 
the hard axis. 
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